Lobe pump system and method of manufacture

ABSTRACT

A method of manufacturing a rotor to be used in a dual-rotor lobe pump system for pumping a material at a periodic rate is provided. The method includes selecting a desired periodic flow rate for the material, selecting a number of lobes for the rotor, and selecting either a thickness of the rotor or a spacing between the dual-rotors&#39; axes of rotation in the lobe pump. The method also includes determining the profile for the rotor based on the desired periodic flow rate, so that when the rotor is operated within the dual-rotor lobe pump system, the material can be pumped at substantially the desired periodic flow rate. In another embodiment of the invention, a lobe pump rotor profile is formed by the method described above.

CROSS-REFERENCE TO RELATED APPLICATION

This application claims the benefit of Provisional Application Ser. No. 60/563,436, filed Apr., 19, 2004, entitled FLOWRATE SYNTHESIS OF LOBE PUMPS, the entire disclosure of which is incorporated herein by reference.

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT

The United States government has certain rights in this invention pursuant to Grant No. CMS-9812847, awarded by the National Science Foundation.

BACKGROUND OF THE INVENTION

Positive displacement rotary pumps, known as “lobe pumps,” are widely used in industries such as pulp and paper, chemical, equipment, food, beverage, pharmaceutical, and biotechnology. Lobe pumps can pump a wide variety of materials at continuous or intermittent flows.

A standard three-lobe pump is shown in FIGS. 1A-1C. Two identical rotors 10, 12 rotate in opposite directions around their respective axes of rotation 14, 16 to mesh as shown. The axes of rotation 14, 16 are separated by a distance l.

Each rotor has multiple lobes 20. The lobes of each of the rotors 10, 12 come in close proximity to the other rotor and to the interior of the lobe pump casing 30, so that material 40 can be trapped between the lobes 20 of the rotors 10, 12 and the pump casing 30.

As the rotors rotate within the lobe pump casing 30, material 40 flows into an inlet end 32 of the casing 30 (FIG. 1A), is subsequently trapped between the lobe 20 of a rotor 10 and the casing 30 (FIG. 1B), and then is pushed out of the pump through the outlet end 34 (FIG. 1C). As the lobes rotate, the material 40 travels around the outside of the rotors 10, 12.

The rotors of a standard lobe pump can be rotated by a driving gear 52 and a driven gear 50, as shown in FIG. 2A. As shown, the rotors 10′, 12′ can each have two lobes 20′ instead of the three shown in FIGS. 1A-1C, or rotors can alternatively be designed to have any number of lobes. The rotor frequency n is the same as the frequency of its driving motor, and is related to a pumping period T by the following expression: ${n = \frac{1}{2{NT}}},$ , where N is the number of lobes on each rotor.

Profiles for the rotors within a lobe pump can be designed using the “deviation function method.” See, e.g., Yang, Tong, and Lin, “Deviation-Function Based Pitch Curve Modification for Conjugate Pair Design,” J. of Mech. Des. v. 121, pp. 579-586 (1999), the entire contents of which are incorporated herein by reference. This method uses a function that describes the deviation of the conjugate pair (or rotor pair) from the profile of a pitch pair, such as a pair of ellipses or circles rotating in opposite directions while maintaining contact. This method allows one skilled in the art to generate a profile of a conjugate pair with a desired geometry so that it matches the rotation of a given pitch pair. For example, the deviation function method could generate a rotor profile with a desired number of lobes of a desired length and noncircularity, etc., that rotates with another rotor similarly to a pair of oppositely rotating circles. This reference allows a broad range of rotor profiles to be generated that correspond to given pitch pairs, but suggests no particular geometry for the rotor or the effects of such geometry.

There are typically two types of lobe pumps used in the industry: conventional, involute lobe pumps and epitrochoidal lobe pumps. FIG. 3A shows a profile of a conventional involute lobe pump rotor. Involute lobe pump rotors have a smooth, continuous profile.

Epitrochoidal lobe pumps have rotors with profiles composed of circular arcs and epitrochoidal curves that do not have first order continuity at some intersections of curve segments. An example of lobe profiles of epitrochoidal rotors is shown in FIG. 4.

Resultant flow rates of conventional lobe pump systems or systems with rotor profiles generated through the deviation function method, described above, have also been previously described by Applicants in “The specific flowrate of deviation function based lobe pumps—derivation and analysis,” Mechanism and Machine Theory 37, pp. 1025-1042 (2002), the entire contents of which are incorporated herein by reference.

In this reference, a normalized flow rate can be derived from a given profile that deviates from an non-circular or circular pitch profile according to a given deviation function, e(θ). Specifically, a flow rate in terms of an angle of rotation θ of the rotor can be expressed as: ${{F(\theta)} = \frac{{\overset{.}{\theta} \cdot {l\left( {b^{2} - {r\left( {l - r} \right)} - {e(\theta)}^{2}} \right)}}w}{2\left( {l - r} \right)}},$ where, referring to FIGS. 2B-C, l represents the distance between the rotors' axes of rotation 140, 160, w is the rotor thickness, b is the lobe length, r is the distance from the axis of rotation 160 of the rotor 120 to a contact point P. The contact point P is the point of contact of the rotors' 140, 160 respective pitch profiles p₁, p₂. e(θ) is the deviation function, or a function showing the deviation of the profile of the actual rotor 120 from its corresponding pitch profile p₁.

It is known that a flow rate of material out of a conventional, involute lobe pump will be a periodic, parabolic function of the angular position θ of the pump rotors, as shown in FIG. 3B. See, Mimmi, 1992; Mimmi and Pennacchi, 1994. The amplitude variation of the periodic function is due to the change of the contact point position of the rotors during the meshing. These periodic functions are described in more detail in, e.g., Yang and Tong, 2000; Bidhendi et al., 1983; and Iyoi and Togashi, 1963. It is also known that the flow rate of material out of epitrochoidal lobe pumps is constant. See Mimmi and Pennacchi, 1994.

One problem present with both existing conventional lobe pump systems is that a user is limited to either a specific constant or a specific periodic parabola flow rate, depending on the type of conventional rotor the user chooses. If a particular periodic flow rate is required for an application, such as a volume of flow that varies sinusoidally with time or angle of rotation, neither of the conventional lobe pump types would be sufficient. Further, even if a periodic parabola or constant type flow rate is required, a user is currently limited to a small number of standard lobe profiles from which to choose. Thus, a user would likely need to employ an entirely different, and costlier, type of pump to achieve a desired flow rate.

SUMMARY OF THE INVENTION

The instant invention is directed to a method of generating lobe pump rotor profiles based on a desired resultant flow rate. One embodiment of the invention provides a method of manufacturing a rotor to be used in a dual-rotor lobe pump system for pumping a material at a periodic rate. The method includes selecting a desired periodic flow rate for the material, selecting a number of lobes for the rotor, and selecting either a thickness of the rotor or a spacing between the dual-rotors' axes of rotation in the lobe pump. The method also includes determining the profile for the rotor based on the desired periodic flow rate, so that when the rotor is operated within the dual-rotor lobe pump system, the material may be pumped at substantially the desired periodic flow rate. In another embodiment of the invention, a lobe pump rotor is formed by the method described above.

BRIEF DESCRIPTION OF THE DRAWINGS

These and other features and advantages of the invention will become better understood when considered in conjunction with the following detailed description and by referring to the appended drawings, wherein

FIG. 1A is a cross section of a conventional lobe pump as material is entering the lobe pump chamber;

FIG. 1B is a cross section of the lobe pump shown in FIG. 1A as the material moves through the chamber;

FIG. 1C is a cross section of the lobe pump shown in FIGS. 1A and 1B as the material begins to move out of the chamber;

FIG. 2A is a diagram of a conventional lobe pump system;

FIG. 2B is a detailed plan view of a pair of lobe pump rotor profiles, generated by the deviation function method, that correspond to non-circular pitch profiles;

FIG. 2C is a detailed plan view of a rotor profile shown in FIG. 2B, showing its corresponding non-circular pitch profile and deviation function;

FIG. 3A is a plan view of a conventional involute lobe pump rotor profile;

FIG. 3B is a diagram of the resultant flow rate of a lobe pump having rotors shaped as shown in FIG. 3A;

FIG. 4 is a plan view of a pair of conventional, epitrochoidal lobe pump rotor profiles;

FIG. 5 is a step diagram of a method according to one embodiment of the invention;

FIG. 6 is a step diagram of a method according to another embodiment of the invention;

FIG. 7A is a plan view of a pair of rotor profiles designed to produce a sinusoidal flow rate, according to one embodiment of the invention;

FIG. 7B is a diagram of the resultant flow rate of a lobe pump having the rotor profiles shown in FIG. 7A;

FIG. 8A is a plan view of a pair of rotor profiles designed to produce a fourth order polynomial flow rate, according to another embodiment of the invention;

FIG. 8B is a diagram of the resultant flow rate of a lobe pump having the rotor profiles shown in FIG. 8A;

FIG. 9A is a plan view of a pair of rotor profiles designed to produce a linear flow rate, according to another embodiment of the invention;

FIG. 9B is a diagram of the resultant flow rate of a lobe pump having the rotor profiles shown in FIG. 9A;

FIG. 10 is a plan view of a pair of rotor profiles designed to produce a constant flow rate;

FIG. 11 is a step diagram of another embodiment of a method according to the invention; and

FIG. 12 is a step diagram of yet another embodiment of a method according to the invention.

DETAILED DESCRIPTION

The instant invention is directed to the design and manufacture of lobe pump profiles that will result in a desired flow rate of material. Referring to FIG. 5, a method for designing a profile includes selecting a desired periodic flow rate for the material. A user may have a particular flow rate function that is required for the application, or the user may merely need a particular maximum flow rate, minimum flow rate, function type (such as parabolic, sinusoidal, polynomial, linear, constant, etc.), and period. A flow rate function may be in many different forms, but the flow of material expressed either in terms of time t or the angle of a rotor's rotation θ will be addressed more specifically below.

A number of lobes for the rotor is then selected, along with a thickness of the rotor or a spacing between the dual rotors' axes of rotation in the lobe pump. The profile is then determined based on the desired periodic flow rate. The determination of the profile can be accomplished by reversing the deviation function method to begin with a desired periodic flow rate and ending with a rotor profile that accomplishes that flow rate.

With reference to FIGS. 2B-3A and 6, another embodiment of the method is described. In this embodiment, the desired flow rate is expressed as a maximum flow rate F_(max), a minimum flow rate F_(min), a function type with some unknown variables F(θ), and a period T. The number of lobes on each rotor is selected to be N, and the distance between the two rotors' axes of rotation is selected to be l.

The function F(θ) of the actual, non-normalized desired flow rate in terms of the angle of rotation θ of the rotor 120 is then generated through known methods using boundary conditions of F(0)=F_(min), F(φ)=F_(max), and ${{\frac{\mathbb{d}{F(\theta)}}{\mathbb{d}\theta}❘_{\phi}} = 0},$ where φ is the angle θ where the pitch profile intersects the generated rotor profile. For circular pitch profiles p, such as that shown in FIG. 3A, $\phi = {\frac{\pi}{2N}.}$

With this function F(θ), and the selected F_(max), F_(min), T, l, and N, half of one lobe profile g is designed according to the following two equations: $\begin{matrix} {{g_{x} = {{\frac{l}{2}\cos\quad\theta} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\cos\left( {\theta + {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)} + \pi} \right)}}}}{{g_{y} = {{{\frac{l}{2}\sin\quad\theta} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\sin\left( {\theta + {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)} + \pi} \right)}\quad{for}\quad 0}} \leq \theta \leq \frac{\pi}{2N}}},{{{{where}\quad{F^{\prime}(\theta)}} = \frac{\mathbb{d}{F(\theta)}}{\mathbb{d}\theta}};{and}}}} & \left. 1 \right) \\ {{g_{x} = {{\frac{l}{2}\cos\quad\theta} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\cos\left( {\theta - {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)}} \right)}}}}{g_{y} = {{{\frac{l}{2}\sin\quad\theta} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\sin\left( {\theta - {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)}} \right)}\quad{for}\quad\frac{\pi}{2N}}} \leq \theta \leq {\frac{\pi}{N}.}}}} & \left. 2 \right) \end{matrix}$

The other half of the lobe profile is then designed to be symmetric to the profile generated by the equations above. Identical lobes can then be designed for a total of N lobes per rotor, which are spaced equally from each other and projecting radially from the axis of rotation 160.

A thickness w of the rotor can be determined according to ${wl}^{2} = {\frac{{NTF}_{\min}^{2}}{\pi\left( {F_{\max} - F_{\min}} \right)}.}$ Alternatively, a desired thickness can be selected and the distance l can be determined through this same calculation. The distance l can then be used to calculate the half lobe profile, as above, and the other half of the lobe profile is then designed to be symmetric to the generated profile.

Although this embodiment is based on generation of a rotor profile that corresponds to a circular pitch profile p (FIG. 3A), rotor profiles may alternatively be generated that correspond to non-circular pitch profiles, such as is shown in FIGS. 2B-2C.

One example of a generation of F(θ) from a function type with unknown variables will now be described. In this example, the function type is selected to be sinusoidal, which can be represented by F(θ)=A₀+A cos αθ, where A₀, A, and α are unknown constants and θ is the angle of rotation of the rotor.

In this case, a normalized function of F(θ) becomes ${{f(\theta)} = {\frac{\pi}{{\pi\quad h^{2}} + {2h}}\left( {{\left( {h + \frac{3}{2}} \right)\left( {h - \frac{1}{2}} \right)} - {\left( {h - \frac{1}{2}} \right)^{2}\cos\quad 2N\quad\theta}} \right)}},{{{where}\quad h} = {\frac{F_{\max}}{F_{\min}} - {{.5}.}}}$ A deviation function can then be determined to be e(θ)=l(h−0.5)cos Nθ. This deviation function can then be inserted into the equation ${F(\theta)} = \frac{{\overset{.}{\theta} \cdot {l\left( {b^{2} - {r\left( {l - r} \right)} - {e(\theta)}^{2}} \right)}}w}{2\left( {l - r} \right)}$ as taught in the prior art and simplified for a circular pitch profile, where $r = {\frac{l}{2}.}$ Further, the function F(θ) can be put in terms of l by substituting the thickness w according to the relation ${wl}^{2} = {\frac{{NTF}_{\min}^{2}}{\pi\left( {F_{\max} - F_{\min}} \right)}.}$ F(θ) is then calculated to be F(θ)=F_(max)−(F_(max)−F_(min))cos² Nθ, which is in the form F(θ)=A₀+A cos αθ through the relation, ${\cos^{2}N\quad\theta} = {\frac{1}{2}{\left( {{\cos\quad 2N\quad\theta} + 1} \right).}}$ If N=2 lobes are selected, the resultant lobe profiles for the desired sinusoidal flow rate type are shown in FIG. 7A. The resultant flow rate in terms of angular position of the rotor is shown in FIG. 7B. As shown, the flow rate varies in amplitude according to the ratio of F_(max) to F_(min), or h+0.5.

Although a sinusoidal function type is discussed above, the function type can alternatively be selected as polynomial, linear, constant, parabolic, and any other continuous functions, and represented as a corresponding function F(θ). Examples of polynomial, linear, and constant flow profiles and their corresponding flow rates in terms of angular rotation of the rotor are shown in FIGS. 8A-B, 9A-B, and 10, respectively.

With reference to FIG. 11, a second embodiment of the method is described. In this embodiment, a desired flow rate is expressed as a function of time, F(t). A number of lobes N and the distance between the axes of rotation l are selected as above.

In this embodiment, F_(max), F_(min), and period T are calculated through known methods from F(t), and a half lobe profile g is designed according to the following two equations: $\begin{matrix} {{g_{x} = {{\frac{l}{2}\cos\quad\frac{\pi \cdot t}{TN}} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(t)}} \right)}{\cos\left( {\frac{\pi \cdot t}{TN} + {\sin^{- 1}\left( \frac{{- {{TNF}^{\prime}(t)}}\sqrt{F_{\max} - F_{\min}}}{{\pi \cdot F_{\min}}\sqrt{F_{\max} - {F(t)}}} \right)} + \pi} \right)}}}}{g_{y} = {{\frac{l}{2}\sin\frac{\pi \cdot t}{TN}} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(t)}} \right)}{\sin\left( {\frac{\pi \cdot t}{TN} + {\sin^{- 1}\left( \frac{{- {{TNF}^{\prime}(t)}}\sqrt{F_{\max} - F_{\min}}}{{\pi \cdot F_{\min}}\sqrt{F_{\max} - {F(t)}}} \right)} + \pi} \right)}}}}{{{{for}\quad 0} \leq t \leq \frac{T}{2}},{{{{where}\quad{F^{\prime}(t)}} = \frac{\mathbb{d}{F(t)}}{\mathbb{d}t}};{and}}}} & \left. 1 \right) \\ {{{g_{x} = {{\frac{l}{2}\cos\quad\frac{\pi \cdot t}{TN}} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(t)}} \right)}{\cos\left( {\frac{\pi \cdot t}{TN} + {\sin^{- 1}\left( \frac{{- {{TNF}^{\prime}(t)}}\sqrt{F_{\max} - F_{\min}}}{{\pi \cdot F_{\min}}\sqrt{F_{\max} - {F(t)}}} \right)}} \right)}}}}{g_{y} = {{\frac{l}{2}\sin\frac{\pi \cdot t}{TN}} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(t)}} \right)}{\sin\left( {\frac{\pi \cdot t}{TN} + {\sin^{- 1}\left( \frac{{- {{TNF}^{\prime}(t)}}\sqrt{F_{\max} - F_{\min}}}{{\pi \cdot F_{\min}}\sqrt{F_{\max} - {F(t)}}} \right)}} \right)}}}}{{{for}\quad\frac{T}{2}} \leq t \leq {T.}}}\quad} & \left. 2 \right) \end{matrix}$

The profile of the other half of the lobe, the remaining lobes, and the rotor thickness are then designed as described above.

Another embodiment of the method is shown in FIG. 12. In this embodiment, a desired flow rate is expressed as a function of the angle of rotor rotation, F(θ). The number of lobes N and distance between the axes of rotation l, is selected as above. F_(max), F_(min), and period T are calculated through known methods from F(θ), and a half lobe profile g is designed according to the following two equations: $\begin{matrix} {{g_{x} = {{\frac{l}{2}\cos\quad\theta}\quad + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\cos\left( {\theta + {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)} + \pi} \right)}}}}{g_{y} = {{\frac{l}{2}\sin\quad\theta} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\sin\left( {\theta + {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)} + \pi} \right)}}}}{{{{for}\quad 0} \leq \theta \leq \frac{\pi}{2N}},{{{{where}\quad{F^{\prime}(\theta)}} = \frac{\mathbb{d}{F(\theta)}}{\mathbb{d}\theta}};{and}}}} & \left. 1 \right) \\ {{{g_{x} = {{\frac{l}{2}\cos\quad\theta} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\cos\left( {\theta - {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)}} \right)}}}}{g_{y} = {{\frac{l}{2}\sin\quad\theta} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\sin\left( {\theta - {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)}} \right)}}}}{{{for}\quad\frac{\pi}{2N}} \leq \theta \leq {\frac{\pi}{N}.}}}\quad} & \left. 2 \right) \end{matrix}$

The profile of the other half of the lobe, the remaining lobes, and the rotor thickness are then designed as described above.

In one embodiment, after the rotor profiles are determined, two identical rotors are formed through conventional methods with a thickness w. The rotors are then placed in a lobe pump on parallel axes of rotation at a distance l from each other. The rotors are then driven by conventional means at a frequency of n=½ NT, where N is the number of lobes and T is the period.

The invention has been described and illustrated by exemplary and preferred embodiments, but is not limited thereto. Persons skilled in the art will appreciate that a number of modifications can be made without departing from the scope of the invention, which is limited only by the appended claims and equivalents thereof. 

1. A method of manufacturing a rotor to be used in a dual-rotor lobe pump system for pumping a material at a desired periodic rate, the method comprising: selecting a desired periodic flow rate for said material; selecting a number of lobes for the rotor; selecting one of a thickness of the rotor or a spacing between the axis of rotation of the rotor and the axis of rotation of an adjacent rotor in said lobe pump; and determining a profile for the rotor based on the desired periodic flow rate, the number of lobes, and the one of the thickness or the spacing, such that when the rotor is operated within said dual-rotor lobe pump system, said material can be pumped at substantially the desired periodic flow rate.
 2. The method of claim 1, wherein the determining a profile further comprises determining a maximum flow rate, a minimum flow rate, and a period from the desired periodic flow rate, the method further comprising: calculating the other of the thickness or the spacing of the rotor based on the one of the thickness or the spacing of the rotor, the maximum flow rate, the minimum flow rate, the number of lobes, and the period; and forming the rotor having the profile, the number of lobes, and the thickness.
 3. The method of claim 2, further comprising forming a second rotor having the profile, the number of lobes, and the thickness.
 4. The method of claim 1, wherein the selecting the desired periodic flow rate comprises selecting a flow rate function F(t) in terms of time t.
 5. The method of claim 4, wherein the determining comprises: determining a maximum flow rate F_(max), a minimum flow rate F_(min), and a period T from the desired periodic flow rate; designing the other of the thickness w or the spacing l of the rotor to satisfy the following formula: ${{wl}^{2} = \frac{{NTF}_{\min}^{2}}{\pi\left( {F_{\max} - F_{\min}} \right)}},$  where N is the number of lobes on the rotor; designing a profile, g_(x) and g_(y), of a first half of one lobe of the rotor to satisfy the following formulas: $g_{x} = {{\frac{l}{2}\cos\quad\frac{\pi \cdot t}{TN}} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(t)}} \right)}{\cos\left( {\frac{\pi \cdot t}{TN} + {\sin^{- 1}\left( \frac{{- {{TNF}^{\prime}(t)}}\sqrt{F_{\max} - F_{\min}}}{{\pi \cdot F_{\min}}\sqrt{F_{\max} - {F(t)}}} \right)} + \pi} \right)}}}$ $g_{y} = {{\frac{l}{2}\sin\frac{\pi \cdot t}{TN}} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(t)}} \right)}{\sin\left( {\frac{\pi \cdot t}{TN} + {\sin^{- 1}\left( \frac{{- {{TNF}^{\prime}(t)}}\sqrt{F_{\max} - F_{\min}}}{{\pi \cdot F_{\min}}\sqrt{F_{\max} - {F(t)}}} \right)} + \pi} \right)}}}$ ${{{for}\quad 0} \leq t \leq \frac{T}{2}},{{{{where}\quad{F^{\prime}(t)}} = \frac{\mathbb{d}{F(t)}}{\mathbb{d}t}};{and}}$ $g_{x} = {{\frac{l}{2}\cos\quad\frac{\pi \cdot t}{TN}} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(t)}} \right)}{\cos\left( {\frac{\pi \cdot t}{TN} - {\sin^{- 1}\left( \frac{{- {{TNF}^{\prime}(t)}}\sqrt{F_{\max} - F_{\min}}}{{\pi \cdot F_{\min}}\sqrt{F_{\max} - {F(t)}}} \right)}} \right)}}}$ $g_{y} = {{\frac{l}{2}\sin\frac{\pi \cdot t}{TN}} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(t)}} \right)}{\sin\left( {\frac{\pi \cdot t}{TN} - {\sin^{- 1}\left( \frac{{- {{TNF}^{\prime}(t)}}\sqrt{F_{\max} - F_{\min}}}{{\pi \cdot F_{\min}}\sqrt{F_{\max} - {F(t)}}} \right)}} \right)}}}$ ${{{{or}\quad\frac{T}{2}} \leq t \leq T};}\quad$ designing a profile of the second half of the one lobe of the rotor, such that said one lobe is symmetrical; and designing a profile of any remaining lobe or lobes of the rotor such that the remaining lobe or lobes have the same profile as said one lobe and all of the lobes are evenly spaced from each other.
 6. The method of claim 1, wherein the selecting the desired periodic flow rate comprises selecting a flow rate function in terms of angular rotation of the rotor F(θ).
 7. The method of claim 6, wherein the determining comprises: determining a maximum flow rate F_(max), a minimum flow rate F_(min), and a period T from the desired periodic flow rate; designing the other of the thickness w or spacing l of the rotor to satisfy following formula: ${{wl}^{2} = \frac{{NTF}_{\min}^{2}}{\pi\left( {F_{\max} - F_{\min}} \right)}},$ where N is the number of lobes; designing a profile, g_(x) and g_(y), of a first half of one lobe of the rotor to satisfy the following formulas: $g_{x} = {{\frac{l}{2}\cos\quad\theta}\quad + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\cos\left( {\theta + {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)} + \pi} \right)}}}$ $g_{y} = {{\frac{l}{2}\sin\quad\theta} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\sin\left( {\theta + {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)} + \pi} \right)}}}$ ${{{for}\quad 0} \leq \theta \leq \frac{\pi}{2N}},{{{{where}\quad{F^{\prime}(\theta)}} = \frac{\mathbb{d}{F(\theta)}}{\mathbb{d}\theta}};{and}}$ $g_{x} = {{\frac{l}{2}\cos\quad\theta} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\cos\left( {\theta - {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)}} \right)}}}$ $g_{y} = {{\frac{l}{2}\sin\quad\theta} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\sin\left( {\theta - {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)}} \right)}}}$ ${{{{for}\quad\frac{\pi}{2N}} \leq \theta \leq \frac{\pi}{N}};}\quad$ designing a profile of the second half of the one lobe of the rotor, such that the one lobe is symmetrical; and designing a profile of any remaining lobe or lobes of the rotor such that the remaining lobe or lobes have the same profile as the one lobe and all of the lobes are evenly spaced from each other.
 8. The method of claim 1, wherein the selecting the desired periodic flow rate comprises: selecting a maximum flow rate F_(max), and a minimum flow rate F_(min); selecting a function type; and selecting a period T.
 9. The method of claim 8, wherein the determining comprises: designing the other of the thickness w or spacing l of the rotor to satisfy following formula: ${{wl}^{2} = \frac{{NTF}_{\min}^{2}}{\pi\left( {F_{\max} - F_{\min}} \right)}},$ where N is the number of lobes; selecting a flow rate function F(θ) having the selected maximum flow rate, minimum flow rate, period, and function type, where θ is an angle of rotation of the rotor; designing a profile, g_(x) and g_(y), of a first half of one lobe of the rotor to satisfy the following formulas: $g_{1x} = {{\frac{l}{2}\cos\quad\theta}\quad + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\cos\left( {\theta + {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)} + \pi} \right)}}}$ $g_{1y} = {{\frac{l}{2}\sin\quad\theta} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\sin\left( {\theta + {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)} + \pi} \right)}}}$ ${{{for}\quad 0} \leq \theta \leq \frac{\pi}{2N}},{{{{where}\quad{F^{\prime}(\theta)}} = \frac{\mathbb{d}{F(\theta)}}{\mathbb{d}\theta}};{and}}$ $g_{1x} = {{\frac{l}{2}\cos\quad\theta} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\cos\left( {\theta - {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)}} \right)}}}$ $g_{1y} = {{\frac{l}{2}\sin\quad\theta} + {\frac{l}{F_{\min}}\sqrt{\left( {F_{\max} - F_{\min}} \right)\left( {F_{\max} - {F(\theta)}} \right)}{\sin\left( {\theta - {\sin^{- 1}\left( \frac{{- {F^{\prime}(\theta)}}\sqrt{F_{\max} - F_{\min}}}{F_{\min}\sqrt{F_{\max} - {F(\theta)}}} \right)}} \right)}}}$ ${{{{for}\quad\frac{\pi}{2N}} \leq \theta \leq \frac{\pi}{N}};}\quad$ designing a profile of the second half of the one lobe of the rotor, such that the one lobe is symmetrical; and designing a profile of any remaining lobe or lobes of the rotor such that the remaining lobe or lobes have the same profile as the one lobe and all of the lobes are evenly spaced from each other.
 10. The method of claim 8, wherein the function type is a continuous function.
 11. The method of claim 10, wherein the continuous function is selected from the group consisting of polynomial, sinusoidal, linear, parabolic, and constant.
 12. A rotor for use in a dual-rotor lobe pump system, made according to the method recited in any one of claims 1-11. 